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An  Almost  Linear-Time  Algorithm 
for  Graph  Realization 

Robert  E.  Bixby  and  Donald  K.  Wagner 

Abstract 

Given  a  {0,  l}-matrix  M,  the  graph-realization  problem  for  M  is  to  find  a  tree  such  that  the 
columns  of  M  are  incidence  vectors  of  paths  in  T ,  or  to  show  that  no  such  T  exists.  An  algo¬ 
rithm  is  presented  for  this  problem  the  time  complexity  of  which  is  very  nearly  linear  in  the 


number  of  ones  in  M. 
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1.  INTRODUCTION 

The  purpose  of  this  paper  is  to  describe  an  efficient  algorithm  for  graph  realization.1  Our 
motivation  is  the  applicability  of  such  an  algorithm  to  finding  network  structure  in  linear  pro¬ 
grams  (see  [2,11]).  Graph  realization  also  provides  a  generalization  of  planarity  testing. 

Fujishige  [7]  has  independently2  given  an  algorithm  with  essentially  the  same  time  bound  as 
that  given  here.  However,  many  of  the  details  in  [7]  are  left  to  the  reader  making  it  difficult  to 
understand  and  verify  the  algorithm.  The  algorithm  in  [7]  and  the  one  given  here  are  similar  in 
that  they  both  are  based  on  an  exponential-time  procedure  proposed  in  [13].  Fujishige  makes  this 
procedure  polynomial  by  employing  the  PQ-tree  data  structure  of  [4],  The  algorithm  given  here 
does  not  use  any  such  special  data  structure,  using  instead  a  graph-decomposition. 

There  are  several  equivalent  definitions  of  the  graph-realization  problem.  The  simplest  is: 
Given  a  {0,  l}-matrix  M,  determine  if  a  tree  T  exists  such  that  the  columns  of  M  are  incidence 
vectors  of  paths  of  T,  and  if  so  find  such  a  tree.  (Note:  It  is  not  necessary  that  all  paths  of  T 
occur  as  columns  in  M,  only  that  all  columns  of  M  occur  as  paths  in  T .)  If  T  exists,  M  is  called 
graphic. 

The  first  polynomial-time  algorithm  for  graph  realization  was  given  by  Tutte  [19].  Other 
polynomial  algorithms  are  given  in  [2,5,7,9,11,14,16,18,20],  Our  algorithm  and  the  algorithm  in  [7] 
have  the  property  that  for  an  M  that  is  not  “entirely  graphic,”  a  maximal  subset  of  “graphic 
columns”  is  produced.  This  property  is  important  in  the  application  to  linear  programming. 

The  remainder  of  the  paper  is  arranged  as  follows.  Section  2  supplies  several  definitions.  In 
section  3,  Lofgren’s  procedure  is  given  along  with  an  informal  discussion  of  how  to  make  it  poly¬ 
nomial.  In  sections  4  and  5  the  main  subroutines  are  developed  in  detail  --  TYPING,  HYPO- 
PATH  and  UPDATE.  In  section  6  the  main  algorithm  GRAPH  is  stated.  Section  7  contains  the 
derivation  of  the  complexity  of  the  algorithm.  Finally,  an  example  is  given  in  the  Appendix. 

'The  paper  is  based  on  the  Ph.D.  dissertation  [22]  of  the  second  author.  An  outline  of  the  main  algorithm  has  also 
recently  been  given  in  [1], 

Preliminary  results  of  the  current  paper  were  presented  at  the  April  30-May  2,  1979  TIMS/ORSA  national  meet¬ 
ing. 
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2.  DEFINITIONS 

We  assume  a  basic  knowledge  of  graph  theory.  See  [3]  for  an  introduction. 

Let  G  =  (N,E)  be  a  graph  with  node-set  N  and  edge-set  E.  An  edge  with  identical  ends  is 
a  loop,  and  two  non-loop  edges  with  the  same  ends  are  parallel.  Let  A  be  a  set  with  AGE.  The 
edge-induced  subgraph  G\A\  is  the  subgraph  with  edge-set  A  and  node-set  consisting  of  those 
nodes  meeting  A.  G\A  is  the  spanning  subgraph  obtained  by  deleting  the  edges  of  A.  A  path  is 
a  “simple  path”  (one  without  repeated  nodes),  and  a  cycle  is  a  “simple  closed  path.”  We  identify 
paths  and  cycles  with  their  underlying  edge-sets.  A  polygon  is  a  connected  graph  the  edge-set  of 
which  is  a  cycle  with  at  least  three  edges.  A  connected,  loopless  graph  on  two  nodes  is  a  bond. 

Let  T  =  (N,A)  be  a  directed  graph  (digraph)  with  arc-set  A  =A(T).  T  is  an  arborescence 
if  its  underlying  graph  is  a  tree,  every  node  except  one  has  indegree  one,  and  the  one  special  node, 
denoted  root(T),  has  indegree  zero.  Let  u  and  v  be  the  tail  and  head,  respectively,  of  some  arc 
of  T.  Then  u  is  the  parent  of  v  and  v  a  child  of  u.  Note  that  every  node  u,  except  root(T), 
has  a  unique  parent,  denoted  here  by  p(u). 

Let  G  =(N,E )  be  a  connected  graph  with  spanning  tree  T.  Then  for  each  edge  e  EE  -T 
the  graph  G[T(J{e}]  contains  a  unique  cycle,  denoted  C(T,e),  the  fundamental  cycle  of  T 
determined  by  e.  Define  a  {0,  l}-matrix  M  as  follows.  Index  the  columns  of  M  on  E-T,  the 
rows  of  M  on  T  and  define  my*,  the  (;,fc)-entry  of  M,  by 

n  if  jec(T,k) 

mi-k  (0  otherwise  . 

M  is  a  fundamental  matrix  of  G  with  respect  to  T. 

Given  a  {0,  l}-matrix  M,  the  graph-realization  problem  is  to  determine  a  connected  graph 
G ,  if  one  exists,  such  that  M  is  a  fundamental  matrix  for  G .  M  is  then  said  to  be  graphic  and  G 
is  a  realization  of  M.  This  definition  is  equivalent  to  the  one  in  the  Introduction.  For  any 
{0,  l}-matrix  M ,  graphic  or  not,  the  set  of  row  indices  corresponding  to  nonzeros  in  column  k 
together  with  the  column  index  k  will  be  denoted  by  Ck  and  be  referred  to  as  a  “fundamental 


cycle”  of  M. 
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Let  G=(N,E )  be  a  connected  graph,  and  let  {EltE2}  be  a  partition  of  E.  Then,  for 
k  >0,  {El,E2}  is  a  k-separation  of  G  if 

\E,\  >k<  |  E2 1 

|fV(G[£1])nM«[^])l  =k 

For  n  a  positive  integer,  G  is  n-connected,  if  it  has  no  ^-separation  for  k  <  n.  Graphs  having  a 
1-separation  are  separable;  2-connected  graphs  are  nonseparable. 

Let  {Ei,E2}  be  a  2-separation  of  a  nonseparable  graph  G.  Let 
A/- ( G [-Ei] )  P)  N{G[E2))=  {«,  a}.  Define  G'  to  be  the  graph  obtained  by  interchanging  in  G[E x] 
the  incidences  of  the  nodes  u  and  v.  Then  G'  is  said  to  be  obtained  from  G  by  reversing 
In  general,  G "  is  2-isomorphic  to  G  if  G"  is  obtainable  from  G  by  a  sequence  of  sub¬ 
graph  reversals. 

A  parallel  class  of  a  graph  G  is  a  maximal  set  of  edges  B  such  that  G[I?]  is  a  bond.  “The” 
simplification  of  G  is  the  graph  obtained  by  deleting  all  loops  and  all  but  one  edge  in  each  paral¬ 
lel  class.  If  G  is  not  a  polygon  or  bond  and  its  simplification  is  3-connected,  then  it  is  prime. 

The  following  important  theorem  is  due  to  Whitney  [24] . 

(2.1)  Theorem.  Let  G  and  G'  be  nonseparable  graphs  on  the  same  edge-set.  Then  G  and  G' 
are  2-isomorphic  if  and  only  if  they  have  the  same  cycles. 

For  a  short  proof  of  Theorem  2.1,  see  [17]  or  [23], 

3.  AN  OUTLINE  OF  THE  ALGORITHM 

Let  M  be  an  rXc  {0,  l}-matrix.  Assume  M  has  no  row  or  column  consisting  of  zeros  only, 
and  let  R  and  G  be  its  sets  of  row  and  column  indices,  respectively.  Define  a  graph  H  with 
node-set  R  (J  C  and  an  edge  for  each  nonzero  of  M.  The  submatrices  corresponding  to  the  con¬ 
nected  components  of  H  are  called  the  blocks  of  M,  and  M  is  said  to  be  nonseparable  if  it  has  just 
one  block.  The  following  result  is  well  known. 
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(3.1)  Theorem.  M  is  graphic  if  and  only  if  each  block  of  M  is  graphic;  moreover,  if  M  is  non- 
separable  and  graphic  with  realization  G ,  then  G  is  nonseparable. 

Linear-time  algorithms  to  find  the  blocks  of  M ,  linear  in  the  number  of  nonzeros,  are  easy 
to  describe  (e.g.,  see  [2]). 

Let  Mk  denote  the  matrix  made  up  of  the  first  k  columns  of  M,  where  rows  consisting 
entirely  of  zeros  have  been  deleted.  M  is  called  totally  nonseparable  if  Mk  is  nonseparable  for 
1  <  A-  <  c .  For  any  nonseparable  matrix  it  is  easy  to  compute  a  permutation  of  the  columns  such 
that  the  resulting  matrix  is  totally  nonseparable.  For  example,  virtually  any  algorithm  for  com¬ 
puting  the  components  of  H ,  above,  will  so  order  the  columns. 

In  view  of  the  above  remarks  we  assume  henceforth  that  M  is  totally  nonseparable.  We 
also  assume,  for  purely  technical  reasons,  that  the  first  column  of  M  is  a  singleton  (see  the  first 
paragraph  of  section  4). 

Where  Ck  is  the  fundamental  cycle  of  column  k  of  M,  define  Pk  =  Ck  p)  I  U  Cy  >  •  A  set  of 

\j<k  \ 

edges  P  of  a  graph  G  is  a  hypopath  of  G  if  P  is  a  path  in  some  graph  2-isomorphic  to  G .  The 
following  statement,  provable  directly  from  Theorem  2.1,  is  Lofgren’s  “subrearrangement 
theorem”  [13]. 

(3.2)  Theorem.  Let  M  be  a  totally  nonseparable  {0,  l}-matrix  with  c  columns.  Assume  for 
some  1  <  k  <  c  that  Mk  is  graphic  with  realization  Gk.  Then  Mk+i  is  graphic  if  and  only  if  Pk+\ 
is  a  hypopath  of  Gk . 

Based  on  Theorem  3.2,  Lofgren  suggested  the  following  procedure  for  testing  whether  M  is 
graphic.  Clearly  M1  is  graphic.  Suppose  there  exists  a  graph  Gk  that  realizes  Mk .  Further,  sup¬ 
pose  Pk+\  is  a  hypopath  of  Gk .  Then  there  exists  a  graph  Gk  ,  2-isomorphic  to  Gk  such  that  Pk+i 
is  a  path  in  Gk  .  Add  the  edges  of  Gk+1-  Pk+l  to  Gk  so  that  they  form  a  path  between  the  ends 
°f  Pk+i  but  are  not  incident  to  any  other  nodes  of  Gk  .  It  is  straightforward  to  verify  that  the 
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resulting  graph  Gk+1  is  a  realization  of  Mk+l.  If  the  above  procedure  breaks  down  at  any  point,  it 
follows  that  M  is  not  graphic. 

To  implement  Lofgren’s  idea  requires  a  polynomial-time  method  for  constructing  Gk'  from 
Gk .  A  natural  approach  is  to  invoke  some  representation  of  Gk  that  “displays”  all  graphs  2- 
isomorphic  to  Gk.  For  this  representation  we  use  a  graph  decomposition,  called  here  t- 
decomposition.  This  notion  of  decomposition  is  a  small  variation  on  a  theory  developed  by  Tutte 
[21,  Chapter  11].  The  data  structures  needed  to  represent  a  t-decomposition  are  elementary. 

The  remainder  of  the  section  is  devoted  to  describing  an  example  decomposition  and  illus¬ 
trating  pictorially  how  it  can  be  used  to  implement  the  above  procedure. 

A  careful  definition  of  t-decomposition  will  be  given  in  the  next  section.  For  now  we  define 
a  t-decomposition  as  an  arborescence  whose  nodes  are  graphs.  Two  nodes,  that  is  graphs,  are 
adjacent  if  they  have  an  edge  in  common,  called  a  “marker”  edge.  It  is  assumed  that  each  graph 
is  a  bond,  polygon  or  prime.  An  example  t-decomposition  is  shown  in  Figure  la.  The  arrows 
between  graphs  indicate  the  pairings  of  the  marker  edges.  The  root  is  Gx.  The  marker  edges 
themselves  are  indicated  by  oriented  dashed  lines. 


^4 


Figure  la:  A  t-decomposition 
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We  associate  with  each  t-decomposition  a  “merged”  graph  obtained  by  identifying  the  com¬ 
mon  marker  edges  and  then  erasing  them.  The  merged  graph  for  the  graph  in  Figure  la  graph  is 
shown  in  Figure  lb.  The  orientations  on  the  marker  edges  specify  the  orientations  in  which  they 
are  identified. 


Figure  lb:  The  merged  graph  for  Figure  la 

Let  G  be  the  graph  in  Figure  lb  and  let  P  =  {a,b ,c ,d,e ,f  ,g}.  P  plays  the  role  of  PA+1  in 
the  Lofgren  procedure.  Each  edge  of  P  lies  in  one  of  the  graphs  of  the  t-decomposition  D  in  Fig¬ 
ure  la.  Consider  the  t-decomposition  D'  in  Figure  2a. 


Figure  2a:  Altered  t-decomposition 
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D'  is  obtained  from  D  by  “reorienting”  Gx  with  respect  to  G2  (that  is,  by  reorientation  of 
its  marker  edge),  “relinking”  the  polygon  G5  (that  is,  by  reordering  its  edges)  and  reorienting  G9 
with  respect  to  Gg.  The  result  is  that  the  merged  graph  G'  of  D'  is  as  given  in  Figure  2b  and 
has  the  property  that  P  is  a  path.  Thus,  P  is  a  hypopath  of  G .  An  algorithm  HYPOPATH  for 
carrying  out  the  above  kinds  of  operations  --  marker  reorientations  and  polygon  relinkings  --  is 
the  subject  of  section  4. 


Figure  2b:  The  merged  graph  for  Figure  2a 

The  final  step  in  an  efficient  implementation  using  t-decompositions  is  to  give  an  update 
procedure.  (A  polynomial-time  but  inefficient  alternative  would  be  to  recompute  the  decomposi¬ 
tion  from  scratch  [10].)  Thus,  suppose  h  is  an  edge  not  in  G' ,  and  that  h  is  added  to  G'  forming 
a  cycle  with  P .  We  must  construct  the  t-decomposition  for  the  new  graph.  It  is  given  in  Figure 
2c.  It  was  obtained  by  merging  G1;  G2,  G3,  G5',  Ge,  Gg  and  G9,  in  the  process  of  which  the 
piece  Gg"  was  “squeezed”  off.  Section  5  describes  the  general  procedure  for  this  update,  the  algo¬ 
rithm  UPDATE. 


Figure  2c:  Updated  t-decomposition 
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4.  THE  HYPOPATH  PROBLEM 

We  begin  by  describing  a  graph  decomposition.  Let  D  be  a  finite  collection  of  nonseparable 
graphs,  and  let  T  be  a  digraph  having  node-set  D  and  two  nodes  joined  by  an  arc  if  these  nodes 
have  some  common  edge.  The  collection  D  is  a  decomposition  if  T  is  an  arborescence  and  any 
two  members  of  D  have  at  most  one  edge  and  no  nodes  in  common.  If  H,KeD  and 
{e }  =  E(H)  E(K),  then  e  is  a  marker  edge  of  H  and  K.  If  K  —  p(H)  (recall  p(-)  denotes 
“parent  of”),  then  e  is  a  child  marker  of  K  and  the  parent  marker  of  H,  in  this  latter  capacity 
denoted  pm(H).  We  assume  that  root{D)  =  root(T)  is  chosen  to  have  at  least  one  edge  not  in 
any  other  node.  Picking  one  such  edge  arbitrarily,  say  m,  we  define  pm{root{D))  =  m .  In  our 
algorithm  we  always  chose  as  m  the  edge  corresponding  to  the  first,  by  assumption  singleton, 
column  of  M.  This  choice  guarantees  that  m  is  never  in  “P”  (see  the  definition  of  the  hypopath 
problem)  and  that  root(D)  is  always  a  bond  given  that  D  is  a  “(-decomposition,”  and  thus 
simplifies  the  statement  of  (R5)  in  section  5. 

Let  e  —pm(H)  have  ends  ut  and  u2  in  Hi  and  and  v2  in  K  =  p(H)  (so,  H  ^  root(D)). 
An  orientation  of  H  with  respect  to  A  is  a  bijection  /  :{wi,  u2}  —+  { v x ,  v2}.  Given  /,  H  is  said  to 
be  oriented  with  respect  to  K.  To  reorient  H  is  to  change  /.  The  decomposition  D  is  oriented  if 
every  child  member  is  oriented  with  respect  to  its  parent.  The  collection  of  these  orientations  is 
referred  to  as  the  orientation  of  D . 

An  oriented  decomposition  D  is  a  t -decomposition  if  (i)  every  member  of  D  has  at  least 
three  edges  and  is  a  polygon,  bond  or  prime;  (ii)  only  bond  members  of  D  have  edges  parallel  to 
their  parent  marker;  and  (iii)  no  two  polygons  or  bonds  have  a  marker  edge  in  common. 

To  relink  a  polygon  member  of  D  is  to  reorder  its  edges,  thus  producing  a  2-isomorphic 

copy. 

Corresponding  to  the  above  decomposition  there  is  a  composition.  Let  D  be  an  oriented 
decomposition  with  H ,K  ED,  K  —p(H),  and  e  =pm(H).  Let  e  have  ends  u j  and  u2  in  H.  To 
merge  H  with  K  is  to  delete  e  from  both  graphs  and  identify  «,•  with  /(«,)  (»  =1,2),  where  /  is 
the  orientation  of  H  with  respect  to  K.  The  resulting  graph  is  denoted  mj{H,K),  or  simply 


10 


m(H,K).  Note  that  m(H,K)  is  uniquely  determined  up  to  the  names  of  the  nodes  being 
identified. 

Replacing,  in  D,  the  graphs  H  and  K  with  m(H,K)  yields  a  new  decomposition  D'  having 
one  fewer  member,  where  the  orientation  of  D  induces  an  orientation  of  D' .  Repeating  this  pro¬ 
cess  an  additional  \D  |  -2  times  yields  a  single  graph  denoted  m(D).  It  is  straightforward  to  ver¬ 
ify  that  m(D )  does  not  depend  on  the  order  in  which  the  members  of  D  are  merged.  The  graph 
m(D )  is  the  merged  graph  of  D  . 

Let  D  — { Q  }  be  an  oriented  decomposition  with  Q  =  p(Hi)  (1  <  »  <  t ).  Let 

be  the  orientations  of  respectively.  Then  Qj  denotes  m(D),  where  / 

refers  to 

(4.1)  Theorem.  Let  G  and  G'  be  2-isomorphic  nonseparable  graphs  and  let  G  =  Qt 
Assume  that  Q  has  at  least  3  edges,  is  a  polygon,  bond,  or  prime,  and  that  if  Q  is  a  polygon, 
then  for  each  H;,  //,  \{ pm (//,-)}  is  nonseparable.  Then  there  exist  graphs  Q '  ,H{  and  an 

orientation  /'  such  that  Q  is  2-isomorphic  to  Q' ,  H;  is  2-isomorphic  to  H{  (l  <  i  <  t)  and 
G'  =  Q'  i  [H{  Further,  if  Q  is  not  a  polygon,  then  Q'  =  Q  ■ 

Proof.  This  may  be  deduced  from  Theorem  2.1  and  results  in  [6].  In  particular,  we  may  replace 
Q  and  each  of  the  H{  (1  <  i  <  t)  by  its  “standard”  decomposition  in  the  sense  of  [6],  Now,  the 
condition  that  when  Q  is  a  polygon,  each  Hj\{pm{H{)}  is  nonseparable,  implies  that  these 
decompositions  taken  together  form  a  “standard”  decomposition  D  of  G  except  for  some  possible 
adjacent  bonds.  But  it  is  easy  to  see  that  merging  these  bonds  does  not  affect  the  class  of  graphs 
obtainable  as  the  merged  graph  by  first  reorienting  D  and  then  merging.  On  the  other  hand, 
because  G  and  G'  have  the  same  cycles,  it  follows  from  the  uniqueness  result  Theorem  18  of  [6] 
that  merging  the  adjacent  bonds  in  D  does  yield  a  “standard”  decomposition  of  G' .  The 
theorem  now  follows  by  applying  Theorem  2.1  and  using  the  easy  fact  that  any  2-isomorphism  of 
a  graph  can  be  carried  out  as  a  2-isomorphism  on  the  simplification  followed  by  the  reinclusion  of 


the  deleted  edges. 
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Let  P  C.E{m{D)).  We  consider  the  following  hypopath  problem:  Given  a  t-decomposition  D 
and  a  set  P  C  E(m(D)),  find  a  t-decomposition  D' ,  if  one  exists,  such  that  m(D')  is  2-isomorphic 
to  m(D)  and  P  is  a  path  of  m(D'). 

Define  D,  the  reduced  t-decomposition  of  D  with  respect  to  P ,  to  be  the  minimal  decompo¬ 
sition  contained  in  D  and  containing  all  members  of  D  that  have  an  edge  of  P .  Clearly,  P  is  a 
hypopath  of  m(D)  if  and  only  if  P  is  a  hypopath  of  m(D).  Note  that  the  arborescence  structure 
of  the  original  t-decomposition  D  induces  an  arborescence  structure  on  D ,  but  it  need  not  be  the 
case  that  root(D)=  root(D). 

The  following  classification  scheme  is  crucial  to  our  method.  Let  H  be  a  graph,  m  a  dis¬ 
tinguished  edge  (m  is  typically  a  parent  marker)  and  4>=^X  C  E(H)-{m}.  Consider  the  following 
mutually  exclusive  arrangements  for  ( H,X,m ): 

(1)  X(J{m}  is  a  cycle; 

(2)  X  is  a  path  with  m  incident  to  one  end-node  and  one  internal  node; 

(3)  X\J{m)  is  a  path  with  m  an  end-edge; 

(4)  {m}  is  a  path  with  m  not  an  end-edge. 

If  ( H,X,m )  satisfies  (i),  set  A{H,X,m)  —  i ;  otherwise  set  A{H ,X,m)  —  5.  Define  the  type  of 
( H,X,m )  by  T(H ,X,m)  =  min{A(H'  ,X,m):H'  is  2-isomorphic  to  H}.  (H,X,m)  is  good  if 

A(H,X,m)=  T(H,X,m)  <  5. 

(4.2)  Lemma.  If  T(H,X,m)  =  i,  then  A{H,X,m)>i.  If  T{H,X,m)= 1,  then  A{H,X,m)  =  1.  If 
T{H,X,m)—  2,  then  A(H,X,m)  =  2  or  5. 

Proof.  The  result  follows  easily  from  the  definition  and  the  easy  part  of  Theorem  2.1. 

For  the  remainder  of  this  section  the  following  notation  is  used.  D  is  the  reduced  t- 
decomposition  of  D  with  respect  to  P,  and  Q£D  is  fixed.  Each  child  of  Q  is  the  root  of  a 
unique  maximal  t-decomposition  contained  in  D.  Let  Si,...,St  be  these  t-decompositions  and  let 
H{  —  m(Si){l<i<t).  are  the  complete  children  of  Q.  Let  m  —pm{Q)  and  define 
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H=  Qf[Hi,...,Ht\  where  /  is  the  orientation  induced  by  D.  Define  mi  =  pm(Hi)  (1  <  i  <  t), 
X=P  Pi  E(H)  and  Wq  =  (P  P)  E(Q))  (J  (mj,...,  mt}.  Note  that  if  Q  is  a  polygon,  then  } 

(1  <  «  <  t)  is  nonseparable  because  root(S{)  cannot  be  a  polygon.  Hence,  Theorem  4.1  applies. 

Suppose  F  is  a  (-decomposition  contained  in  D  with  K=m(F)  and  P  In 

what  follows  T(K)  T(K,  Pf~]E(K),  pm(root(F))),  A(K)  :=  A (K,  P(~)E(K),  pm(root(F))) 
and  AT  good  means  T{K)  =  A(K)  <  5. 

The  following  three  lemmas  provide  some  basic  results.  The  proofs  are  straightforward  and 
are  left  to  the  reader.  An  end-node  of  H[X ]  is  a  node  incident  to  exactly  one  edge  of  X. 

(4.3)  Lemma.  Suppose  X  is  a  disjoint  union  of  paths  of  H. 

(a)  If  T(H{)  >  1  (respectively,  >  3)  (some  t),  then  contains  an  end-node  (respectively, 
two  end-nodes)  of  i/[X]  not  incident  to  m,  and  T(H)>  1  (respectively,  >3). 

(b)  If  H[X\  has  an  end-node  (respectively,  two  end-nodes)  in  H{  (some  j)  and  not  incident 
to  m,-,  then  T(H{)  >  1  (respectively,  >  3). 

By  Lemma  4.3,  T(H)<  4  implies  at  most  two  of  the  (Hit X(, m,)  are  not  type  1.  Further,  if 
for  some  i,  T[H{)  =  4,  then  T(H})=  1  for  jj^i. 

(4.4)  Lemma.  If  A(H)  <  4,  then  A(Hj)  <  4  (1  <  i  <  ().  Moreover,  if  T(H,)  =  5  for  some  t,  then 
P  is  not  a  hypopath  of  m{D). 

(4.5)  Lemma.  Let  H{  be  a  2-isomorphic  copy  of  Hf  for  some  i.  Assume  A(Hj')  =  A(H{).  Then 
there  is  an  orientation  of  H{  such  that  replacing  Hi  by  H{  leaves  A(H)  unchanged. 

Theorems  4.6  and  4.8  are  the  heart  of  the  algorithm  for  solving  the  hypopath  problem. 

(4.6)  Theorem.  Suppose  that  Q  is  not  a  polygon.  If  each  //,■  is  good,  then  either  there  exists  an 
orientation  /'  such  that  Qj-[Hi,...,Ht]  is  good,  or  T(H)  =  5. 
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Proof.  Assume  T(H)  <  4.  Then  by  Theorem  4.1,  since  Q  is  not  a  polygon,  there  are  2- 
isomorphic  copies  H{  of  respectively,  and  an  /'  such  that  A(H')=  T(H)  <  4, 

where  H'  —  Qf[H{ Ht'].  By  Lemma  4.4,  A(H/)  <  4  (1  <  i  <  t),  and  by  Lemma  4.5  we  may 
assume  not  all  H{  are  good.  Suppose  H{  is  not  good.  Then  by  Lemma  4.2,  T(Hl)  =  3  and 
A(H{)  =  4,  and  so  by  Lemma  4.3,  T(Hi)  =  l  (»  >  1).  Now,  it  is  obvious  that  A(H')>  4,  and 
that  replacing  H{  by  H1}  in  either  orientation,  does  not  decrease  A(H'). 

(4.7)  Computational  Remark.  /'  is  straightforward  to  compute.  This  remark  applies  as  well  to 
Theorems  4.8-4.10.  In  particular,  by  Lemma  4.3,  either  exactly  one  H{  is  of  type  4,  and  all  others 
are  of  type  1,  or  at  most  two  have  type  2  or  3  and  all  others  are  type  1.  In  the  first  case  no 
reorientation  is  needed.  In  the  second  case  there  are  only  two  possible  orientations  for  the  non¬ 
type  1  Ht,  and  in  determining  which  of  these  are  good  only  the  degrees  of  the  ends  of  pm(H{ )  in 
Hi[P(~]E(Hi)]  are  needed.  /'  can  thus  be  found  in  time  linear  in  |  Wq  \ . 

A  similar  theorem  holds  when  Q  is  a  polygon,  but  because  polygons  may  have  nontrivial  2- 
isomorphisms  a  slight  modification  is  needed.  Let  b  be  the  number  of  non-type  1  complete  chil¬ 
dren  of  H.  By  Lemma  4.3  it  may  be  assumed  that  b  <  2.  Let  Z  C  Wq  be  the  set  of  child  mark¬ 
ers  of  H  that  correspond  to  the  non-type  1  complete  children.  Assume  that  either 
Z  =  (j),  =  { m x}  or  Z  =  {m1,m2}.  The  appropriate  relinking  procedure  is: 

Procedure  RELINKl(Q) 

Q  must  be  a  polygon,  and  Z  must  satisfy  \Z\  <2.  Relink  Q  so  that  Wq-Z  is  a  path,  one 
end  of  Wq-Z  is  incident  to  m,  the  other  end  is  incident  to  mj  (if  it  exists),  and  m2  (if  it 
exists)  is  incident  to  m. 

• 

(4.8)  Theorem.  Suppose  RELINK1((J)  has  been  applied.  If  each  //,•  is  good,  then  either  there 
exists  an  orientation  /'  such  that  Qj- [Hl,...,Ht\  is  good  or  T(H)  =  5. 


Proof.  There  are  essentially  four  cases.  If  Z  =  <j),  the  result  is  clear. 
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If  Z—{mx)  and  |  Wq  |  <  |E((j))|-l,  then  T(H)>  3.  If  A(HX)  =  2  or  3,  then  it  is  easily 

seen  that  there  exists  an  /'  such  that  A(Qr[H1,...,Ht})  =  3.  If  A{H1)  =  4,  then  Theorem  4.1 

implies  T(H)  =  5  since  T(H)~>  4  by  Lemma  4.3  and  Wq  does  not  have  enough  edges  for 
H[X\ J  {m}]  to  be  connected. 

If  Z  —  {mi}  and  |  Wq  \  =  \  E(Q)  |  -1,  then  T(H)  >2.  If  A{HX)  =  2,  3  or  4,  then  clearly 

}'  exists  such  that  A{Qj[Hx,  ...,Ht})  —  A(HX),  and  so  if  A(Hl)  =  2  we  are  done.  In  the  case 

A(jfiTi)  =  3,  the  easy  part  of  Theorem  2.1  implies  T(H)  >A(Hi)  since  X[J  {m}  does  not  contain 
a  cycle,  and  when  A(HX)  =  4,  Lemma  4.3  implies  T(H )  >A(HX). 

If  Z={m1,  m2}  and  |  Wq  |  <  |E(Q)|-1,  then  T(H)>  4.  On  the  other  hand,  by  Lemma 
4.3,  T(H)  <4  only  if  T(H{)< 3  (»  =  1,2).  In  this  case  T(H)  =  4  is  clearly  achievable. 

Suppose  Z  =  {m1,m2}  and  \WQ\  =  |£(Q)|-1.  Then  T(H)>  4.  If  T(H1)=T{H2)  =  2,  then 
T(H)  =  5  since  X[_J{m}  contains  a  cycle.  If  T(H{)>  4  (*  =  1,2),  then  T(H)  =  5  by  Lemma  4.3. 
In  all  other  cases  it  is  easy  to  see  T(H)  =  4  is  achievable. 


Algorithm  TYPING 

Input:  The  reduced  t-decomposition  D  corresponding  to  a  <-decomposition  D  and  a  nonempty  set 
P  C  E(m(D));  the  depth  partition  *r=(7r0, ..., *r,)  of  D,  where  GEjry(0<  /  if  the  unique  path 
from  G  to  root(D)  has  j  arcs.  We  assume  s  >  1. 

Output:  The  conclusion  that  P  is  not  a  hypopath  of  m{D),  or  a  reoriented  D,  and  hence  a 
reoriented  D,  such  that  each  complete  child  Hi  of  root(D)  is  good  and  T(H{)  is  known. 

Comment.  Steps  Tl  and  T3  are  modified  in  section  5  under  “Rules  for  selecting  KX,K2." 

Step  Tl.  For  each  polygon  HEtt,  apply  RELINKl(tf).  If  T(H)  =  5  for  some  Hen,,  stop  - 
Lemma  4.4  implies  P  is  not  a  hypopath  of  m(D).  (Note  that  each  Hen,  is  a  polygon, 
bond,  or  prime,  and  so  the  computation  of  T(H)  is  easy.  The  data  structures  required  to  do 
this  computation  in  time  linear  in  Pf^\E(H)  are  discussed  in  section  6.)  Set  i*—  «— 1. 

Step  T2.  If  *  =0,  stop  --  the  desired  D  has  been  found. 

Step  T3.  Let  QEnit  let  Hlt...,Ht  be  the  complete  children  of  Q  in  D,  and  let 
H  =  Qf[Hl,...,Ht]  where  /  is  the  current  orientation  of  D.  If  more  than  two  of  the  H  are 
not  type  1,  stop  —  Lemma  4.3  implies  P  is  not  a  hypopath  of  m(D).  If  Q  is  a  polygon, 
apply  RELINKl(<2)-  Find,  if  possible,  an  orientation  /'  and  a  corresponding  new  D  such 
that  Qf[Hi,...,Ht]  is  good  -  see  (4.7).  If  T(H)  =  5,  stop. 

Repeat  the  above  procedure  until  a  stop  occurs,  or  every  Q€*r,-  has  been  treated. 
When  the  latter  occurs,  set  *«—  f-1  and  go  to  Step  T2. 
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The  considerations  required  are  slightly  different  when  Q — root(D)  (although  the  proofs 
are  essentially  the  same).  Typing  is  no  longer  important;  we  wish  rather  to  determine  whether  or 
not  P  can  be  made  into  a  path.  For  example,  in  this  case  T(H)  =  5  is  preferred  to  T(H)  =  4. 

(4.9)  Theorem.  Suppose  Q  =  root(D)  and  Q  is  not  a  polygon.  If  each  H{  is  good,  either  there 
exists  an  orientation  mapping  }'  such  that  P  is  a  path  of  Qy  or  P  is  not  a  hypopath 

of  H. 

For  polygons  the  following  relinking  procedure  is  needed. 

Procedure  RELINEv2(  Q ) 

Q  must  be  a  polygon.  Relink  Q  so  that  Wq-Z  is  a  path  and,  if  m,-  exists  («  — 1,2),  so  that 
m,-  is  incident  to  an  end  of  the  path  Wq-Z. 


(4.10)  Theorem.  Suppose  Q  =  root(D)  and  that  RELINK2(<2)  has  been  applied.  If  each  H{  is 
good,  then  either  there  exists  an  orientation  /'  such  that  P  is  a  path  of  Qj- {Hlt  or  P  is 

not  a  hypopath  of  H. 


Algorithm  HYPOPATH 

Input:  A  t-decomposition  D  and  a  set  P  C  E(m(D)). 

Output:  A  (new)  ^decomposition  D  such  that  P  is  a  path  of  m(D)  and  m(D)  is  2-isomorphic  to 
the  input  m(D),  or  the  conclusion  that  no  such  2-isomorphism  exists,  i.e.,  that  P  is  not  a  hypo¬ 
path  of  D. 

Step  Hi.  Compute  the  reduced  t-decomposition  D  for  P .  If  D  has  just  one  member,  let 
this  member  be  Q  and  go  to  Step  H3. 

Step  H2.  Apply  TYPING.  If  this  results  in  the  conclusion  that  P  is  not  a  hypopath,  stop; 
otherwise,  set  Q  <—  root(D). 

Comment.  Step  H3  is  expanded  in  section  5  under  “Rules  for  selecting 

Step  H3.  Let  be  the  complete  children  of  Q  in  D.  If  more  than  two  H  are  not  of 

type  1,  stop  --  P  is  not  a  hypopath.  Apply  RELINK2(Q)  if  Q  is  a  polygon.  Find,  if  possi¬ 
ble,  an  orientation  f'  and  corresponding  D  such  that  P  is  a  path  of  Qy[Hlt...,  Ht\  —  see 
(4.7). 
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5.  DECOMPOSITION  UPDATE 

Assume  the  following  are  given:  A  f-decomposition  D,  a  path  P  of  m(D),  and  the  reduced 
f-decomposition  D  of  D  with  respect  to  P.  Let  G  be  a  set  such  that  C  (~]E(m{D))  =  P  with 
C-P^4>.  The  purpose  of  this  section  is  to  give  a  method  to  determine  a  ^-decomposition  D *  for 
the  graph  obtained  from  m(D)  by  adding  the  edges  of  C-P  so  that  C  is  a  cycle,  and  C-P  is 
incident  to  m(D)  at  exactly  two  nodes.  (In  terms  of  the  notation  in  section  3,  P  =  P*+1,D  =  ZV, 
G=Gk+ 1  and  D*  =  Dk+1.) 

The  idea  for  finding  D*  is  intuitively  simple.  First  we  select  appropriate  graphs  Kx,K2ED 
containing  the  ends  of  P.  If  it  happens  that  Kl  =  K2,  then  the  construction  of  D  *  changes  D 
very  little.  When  Kl  more  complicated  modifications  are  involved.  Where  R  is  the  unique 

path  in  D  between  K{  and  K 2,  we  first  “squeeze”  off  certain  parts  of  any  polygons  that  occur  in 
R.  Then  the  remaining  graphs  are  merged  and  C-P  appropriately  added. 

We  introduce  a  convention.  For  a  given  node  x  of  m(D),  there  may  be  several  nodes  of 
members  of  D  that  are  merged  into  x.  In  the  remainder  of  this  section  we  identify  all  these  nodes 
with  the  name  x. 

Using  the  following  rules,  the  graphs  Klt  K2  are  conveniently  found  while  executing  TYP¬ 
ING  and  HYPOPATH. 

RULES  FOR  SELECTING  KUK2 

Comment.  The  end-nodes  of  P  contained  in  KltK2  are  computed  while  applying  rules 

(Rl)-(R4).  They  will  be  encountered  naturally  during  the  calculations  in  Steps  Tl,  T3  and 

H3. 

Additions  to  Steps  Tl  &  T8  in  TYPING: 

(Rl)  If  T(H;)  =  1  (all  i)  and  T(H)  =  2  or  3,  either  set  K1-*—Q  if  Kx  is  not  yet  assigned,  or 
set  K2*~  Q  ■ 

(R2)  If  T(H()  =  1  (all  »)  and  T(H)  =  4,  set  KUK2+-Q. 

(R3)  Suppose  T(Hk)  =  2  or  3  for  some  k,  T(Hf)  =  1  (all  i^=k)  and  T(H)  =  4.  Let  K2  be 
the  node  on  the  unique  path  in  D  between  Q  and  Ki  that  is  nearest  Ki  and  contains 
the  same  end  of  P  as  Q . 
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Additions  to  Step  HS  of  HYPOPATH: 

(R4)  If  Ki  has  not  been  assigned,  set  KltK2*—  Q.  Otherwise,  let  K2  be  the  node  on  the 
unique  path  in  D  between  Q  and  Ki  that  is  nearest  K1  and  contains  the  same  end  of 
P  as  Q. 

(R5)  Suppose  KX,K2  have  been  selected  and  KX  =  K2.  If  the  ends  of  P  are  the  ends  of 
pm{Kx)  and  K1  is  not  a  bond,  set  Kx, K2*—  p{Kf).  If  the  ends  of  P  are  the  ends  of 
another  marker  m,-  of  Kx,  Kx  is  a  polygon,  and  m,-  is  an  edge  of  a  bond  G,  set 
KltKa*-G. 

• 

The  application  of  the  final  rule  (R5)  is  needed  for  the  proof  of  Lemma  5.4. 

In  the  above  list  of  rules,  (Rl)  is  justified  by  Lemma  4.3.  In  particular,  (4.3a)  implies  N(H) 
has  an  end-node  x  of  P  not  incident  to  pm(Q),  and  (4.3b)  implies  xEN(Q).  (R2)  is  justified 
similarly.  For  (R3),  first  note  that  Kx  has  already  been  selected  because  (Rl)  must  have  been 
applied  to  Hk  or  one  of  its  descendants.  Now  by  Lemma  4.3,  P  has  an  end  x£N(Q)  not  incident 
to  pm(Q),  which  implies  that  K2  could  not  previously  have  been  assigned.  However,  x  may  be 
incident  to  a  child  marker  of  Q,  in  which  case  K2^Q  is  possible.  Consider  now  (R 4).  Clearly 
N(H)  contains  both  ends  of  P.  If  neither  Kx  nor  K2  has  been  assigned,  then  T(H,)=  1  (all  i). 
Hence,  Lemma  4.3  implies  both  ends  of  P  are  in  N(Q),  as  required.  Similar  reasoning  justifies 
the  choice  of  K2  when  K j  has  already  been  assigned. 

(5.1)  Lemma.  The  application  of  (Rl)  during  the  execution  of  HYPOPATH  and  TYPING  finds 
Kx,K2ED  such  that  «,  eN(iC,)  (*  =  1,2)  where  uuu2  are  the  distinct  end-nodes  of  P  in  m(D). 
Where  R  is  the  unique  path  in  D  between  Kx  and  K2,  no  internal  node  of  R  contains  ux  or  u2. 

Proof.  By  the  remarks  preceeding  the  lemma  we  need  only  prove  the  last  statement.  There  is 
nothing  to  prove  if  Kx  —  K2.  If  (R3)  or  (R4)  was  applied  we  use  the  fact  that  ux  is  not  incident 
to  pm(Ki),  and  so  is  contained  only  in  descendants  of  Kx.  The  only  remaining  possibility  is  that 
(Rl)  has  been  applied  twice.  But  by  (4.3a),  in  this  case  neither  K}-  is  a  descendant  of  the  other, 
and  Uj  (y  =  1,2)  is  contained  only  in  descendants  of  Kj.  This  completes  the  proof. 

Given  that  KltK2  have  been  selected,  we  can  carry  out  the  actual  update.  Recall  that  t- 
decompositions  have  an  underlying  arborescence  structure,  and  are  oriented.  In  what  follows, 
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when  new  members  are  added  to  D  to  form  D *  we  will  not  explicitly  specifiy  the  new  root,  how 
new  arcs  of  the  arborescence  are  directed,  or  how  the  orientations  of  new  members  are  deter¬ 
mined.  The  specification  of  the  new  root  is  straightforward:  It  will  always  be  the  unique  node 
that  contains  the  parent  marker  of  the  old  root.  Given  a  root,  the  directions  of  the  arcs  of  D*  are 
uniquely  specified.  Now,  for  the  orientation  of  D* ,  the  case  KX=K2  is  easy,  and  when  Kxf^K2, 
with  one  exception,  we  employ  the  natural  orientation;  thus,  in  procedure  SQUEEZE  we  choose 
the  orientation  such  that  the  merged  graph  is  the  one  we  started  with.  The  one  exception  occurs 
in  Step  UO  where  the  orientation  of  “{/}  [J  (C-P)”  is  arbitrary,  and  in  any  case  can  only  be 
determined  once  “/”  has  been  added  to  another  member  of  D* . 

We  introduce  the  following  procedure  for  application  in  Step  U2  of  UPDATE.  Note  that 
D  *  and  R  are  defined  in  UPDATE. 

Procedure  SQUEEZE(L) 

L  must  be  a  path  in  some  polygon  S  of  R.  If  L  has  fewer  than  two  edges,  do  nothing. 
Otherwise,  let  /'  be  a  new  element  and  replace  S  in  D'  by  the  two  polygons  formed  by 
adding  /'  to  the  paths  L  and  S-L,  respectively.  Replace  S  in  R  by  the  polygon  formed 
from  S-L. 


Algorithm  UPDATE 

Input:  A  (-decomposition  D,  a  path  P  of  m(D),  and  a  set  C  such  that  C  p|  E(m(D))  =  P  and 
C-Py^<f>;  nodes  KX,K2  of  D  and  nodes  ux,u2  of  KX,K2,  respectively,  such  that  the  conclusions 
of  Lemma  5.1  hold.  We  also  assume  that  (R5)  has  been  applied. 

Output:  A  (-decomposition  D*  of  the  graph  obtained  from  m(D)  by  adding  the  edges  of  C-P  so 
that  C  is  a  cycle,  and  C-P  is  incident  to  m(D)  at  exactly  two  nodes. 

Step  UO.  Set  Z)*<—  D.  If  |  C-P  j  =1  set  {/}•<—  C-P;  otherwise,  let  /  be  a  new  element, 

form  a  polygon  with  edge-set  {/}U(G-P)  (the  order  of  the  edges  is  irrelevant)  and  add 

this  polygon  to  D’ .  If  KX  =  K2,  go  to  Step  Ul;  otherwise,  go  to  Step  U2. 

Step  Ul.  (Case  KX=K2)  Apply  the  appropriate  one  of  (Ul.l)-(U1.3)  and  stop: 

(Ul.l)  If  Kx  is  not  a  polygon,  join  ux  and  u2  by  /  in  Kx. 

(Ul.2)  Suppose  Kx  is  a  polygon  and  u x,u2  are  joined  in  Kx  by  edge  /'.  Let  /"  be  a  new 
element.  Replace  /'  by  f”  in  Kx  and  add  a  bond  with  edge-set  {/,/',/"}  to  D*. 

(Ul.3)  Suppose  Kx  is  a  polygon  and  ux,u2  are  not  adjacent.  Then  there  are  two  distinct 
paths  in  Kx,  say  Lx  and  L2,  joining  ux  and  u2.  Let  fx,f2  be  new  elements.  Delete 
Kx  from  D  *,  add  polygons  formed  by  joining  the  ends  of  L,-  with  /,-(t— 1,2),  and 
add  a  bond  with  edge-set  {y ,  / 1 ,  / 2}  • 
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Step  U2.  (Case  KX^K2)  Let  R  be  the  unique  path  Kx  =  J j,...,  Je+X=K2  in  D *  joining 

K1  and  K2.  Let  {my}  =  F(Jy)  p|  E(JJ+1)  (1  <  j  <  s).  Apply  (U2.1)-(U2.4)  and  stop: 

(U2.1)  Suppose  1  <  j  <  s,  Jj  is  prime  and  {my_!,  my}  is  a  cycle  of  Jy,  or  Jy  is  a  bond  on 
at  least  4  edges  and  p(Jy)  is  not  in  P.  Let  /'  be  a  new  element.  Let  J/  be  Jj 
with  {my_i,  my}  deleted  and  /'  added  (with  the  same  ends),  and  let  B  be  a  bond 
with  edge-set  {my.!,  /',  my}.  Replace  Jj  in  R  by  S;  delete  Jj  from  D*  and 
replace  it  by  Jj'  and  B.  (Note:  There  can  be  at  most  one  Jy  as  above.) 

(U2.2)  If  K j  is  a  polygon,  let  Ll,L2  be  the  two  paths  joining  ml  and  «x.  Apply 
SQUEEZE(L,-)  («  =  1, 2).  Do  the  same  for  K2,  if  K2  is  a  polygon. 

(U2.3)  For  each  internal  Jy  in  R  that  is  a  polygon,  let  Lx,  L2  be  the  two  components  of 
Jj  -  {my,  m;+1}  and  apply  SQUEEZE(L,-)  (»  =1, 2). 

(U2.4)  Let  G  be  m{R)  with  /  joining  ux  and  u2.  Delete  R  from  D'  and  add  G . 


We  must  now  prove  that  D‘  is  a  t -decomposition  and  that  it  is  the  “right”  t- 
decomposition.  The  second  is  the  easier  of  the  two  assertions,  and  is  treated  in  the  following 
lemma. 


(5.2)  Lemma.  After  UPDATE  has  been  applied,  m(D')  equals  the  graph  obtained  from  m(D)  by 
adding  the  edges  C-P  so  that  G  is  a  cycle  and  C-P  is  incident  to  m{D)  at  exactly  two  nodes. 

Proof.  This  proof  reduces  to  a  careful  reading  of  the  steps  of  UPDATE.  We  treat  only  the  case 
|  C-P  |  >1.  In  Step  UO  denote  {/ }  {J  {C -P)  by  S'.  Note  that  S'  is  a  member  of  D* .  Let 
D i*  denote  the  final  D‘  with  S'  deleted  and  /  deleted  from  the  remaining  member  of  D*  that 
contains  it  as  an  edge. 

First,  we  observe  that  m(Di*)  =  m(D).  Step  UO  and  (Ul.l)  present  no  difficulties.  In 
(Ul.2),  once  /  is  deleted,  merging  the  remaining  2-bond  {/',/"}  with  Kx  yields  D;  in  (U1.3), 
merging  the  remaining  2-bond  {/i,/2}  with  the  Lx  and  L2  polygons  yields  D .  Now  consider  Step 
U2.  Deleting  /  in  this  case  just  means  leaving  out  a  part  of  (U2.4).  But  then  it  is  straightfor¬ 
ward  to  see  that  none  of  (U2.1)-(U2.4)  change  m(D*).  This  proves  m{Dx>)  =  m{D). 

Now  to  complete  the  proof,  note  that  /  is  added  in  either  (U1.1)-(U1.3)  or  (U2.4).  It  is  clear 
in  (Ul.l)  and  (U2.4)  that  /  joins  ux  and  v2  in  m(D)=  m(P1‘).  Steps  (Ul.2)  and  (Ul.3)  are  han¬ 
dled  as  in  the  previous  paragraph.  Now,  given  that  /  has  been  added,  it  only  remains  to  reinsert 
S' .  But  clearly  this  has  the  desired  effect  of  adding  C-P  as  a  path  meeting  m(D )  at  ux,  u2. 
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The  main  step  in  proving  that  D*  is  a  (-decomposition  is  the  following  lemma. 

(5.3)  Lemma.  G  in  (U2.4)  is  prime. 

Proof.  Let  x  and  y  be  distinct  nodes  of  G.  We  show  that  there  are  three  node  disjoint  (except 
for  x  and  y)  paths  joining  x  and  y.  Assume  that  in  (U2.4),  xEN(J{)  and  yEN(Jk)  where  %  <  k. 
We  also  assume  {x ,y}  f~)  {«i,«2}  =4>\  the  proof  when  this  latter  assumption  fails  is  an  easy  varia¬ 
tion  of  the  following  arguments. 

By  (U2.2)  and  (U2.3)  we  may  assume  neither  J{  nor  Jk  is  a  polygon.  If  i  =  k,  then  since  Jk 
is  either  prime  or  a  bond  with  at  least  three  edges,  it  contains  three  node-disjoint  paths  from  x  to 
y.  Suppose  one  of  these  paths  L  contains  m,-.  Then  the  merge  of  Jf+i,...,  Jt)  since  this  graph  is 
2-connected,  contains  a  cycle  C  containing  m,-,  and  (G  (J  L)- {m,}  is  a  path  in  G .  Applying  this 
procedure  at  most  twice,  we  obtain  the  desired  paths  in  G. 

Suppose  that  ij^k.  Consider  J{.  If  i  7^  1  and  /,■  is  prime,  has  at  least  three 

nodes  by  (U2.1)  and  there  are  three  node-disjoint  paths  from  x  to  three  of  these  nodes,  including 
the  ends  of  m,-,  and  such  that  mf  and  m,_i  are  not  in  these  paths.  Let  a!  be  the  end  of  m,_!  used 
and  let  a2  be  the  other  end  of  m,-^.  If  *  =  1,  find  three  node-disjoint  paths  in  Jx  from  x  to  ux  and 
the  ends  of  m1  (and  not  containing  mj);  in  this  case  a2  is  undefined  and  ai  =  Ui.  In  the  case  that 
Jj  is  a  bond,  let  o/l  =  x  and  a2  be  the  other  node  in  J; .  Perform  the  symmetric  construction  for 
Jk  defining  j3x  and  /?2.  Let  Lx  be  a  path  in  the  merge  of  from  to  alt  avoiding  a2.  L2 

is  defined  similarly  for  u2,/32.  Finally,  let  G  be  a  cycle  of  the  merge  of  Jk_x  (empty  if 

*  +  !  —  &)  containing  m,-  and  mk_ j.  Then  LIUL21>J  {/}  U(C'-{m,-,mt_1})  together  with  the 
paths  constructed  in  J,-  and  Jk  gives  the  desired  three  paths  in  G . 

We  are  now  prepared  to  prove 

(5.4)  Lemma.  D*  is  a  (-decomposition. 

Proof.  We  must  verify  that  D "  has  the  following  properties:  (a)  It  is  a  decomposition,  (b)  every 
member  has  at  least  three  edges,  (c)  there  are  no  adjacent  polygons  or  bonds,  (d)  no  non-bond  has 
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an  edge  parallel  to  its  parent  marker,  and  (e)  every  member  is  a  bond,  polygon,  or  prime. 

Most  of  the  work  in  verifying  (a)-(e)  involves  simply  a  careful  reading  of  UPDATE,  as  in  the 
proof  of  Lemma  5.2.  We  mention  only  the  highlights. 

For  (a)  note  that  in  Step  UO  and  in  the  case  \G—P\  >1,  the  conditions  for  D‘  to  be  a 
decomposition  are  violated  since  the  polygon  (C-P)\^J  {f}  has  no  edges  in  common  with  any 
other  member  of  D*.  However,  it  is  easy  to  verify  that  /  is  later  added  to  exactly  one  other 
member  of  D  * ,  thus  restoring  the  required  property. 

In  proving  (b)  note  that  \C-P\  >1  implies  \{C-P)  l _ J  {/ }  |  >3  in  Step  UO.  Also  note 

that  whenever  SQUEEZE(L)  is  applied,  no  change  occurs  unless  L  has  at  least  two  edges  ( S-L 
always  has  at  least  two  edges).  Finally,  note  that  since  J,-  in  (U2.1)  is  prime,  or  a  bond  on  at 
least  4  edges,  |  E(Jj)  |  - 1  >  3. 

Part  (e)  is  easy  to  verify  in  view  of  Lemmna  5.3.  Consider  now  (c).  Step  UO  creates  no  new 
adjacencies,  and  so  no  violations  to  (c).  However,  a  polygon  may  be  created  and  it  must  then  be 
verified  that  edge  /  is  added  to  no  other  polygon.  That  this  is  true  follows  from  inspection  and 
Lemma  5.3.  Note  also  that  no  violation  of  (c)  is  created  in  (U1.2).  Here  /'  cannot  be  in  another 
bond  of  D*  by  (R5). 

Finally,  consider  (d).  (Ul.l)  preserves  (d)  by  (R5).  Steps  (Ul.2),  (Ul.3),  and  (U2.1)-(U2.3) 
increase  the  “edge  multiplicities”  only  in  bonds.  In  (U2.4),  «]  and  u2  are  not  adjacent  until  /  is 
added,  and  otherwise  edge  multiplicities  are  not  increased.  Hence,  if  in  (U2.4)  the  parent  marker 
of  G  is  parallel  to  another  edge,  then  it  was  the  parent  marker  of  a  bond  in  R  with  at  least  4 
edges.  This  possibility  is  ruled  out  by  (U2.1).  This  completes  the  proof. 


6.  THE  MAIN  ALGORITHM 

We  begin  by  stating  GRAPH.  This  is  followed  by  a  listing  of  the  data  structures  that  are 
needed  to  implement  the  algorithm,  and  a  discussion  of  how  these  data  structures  are  updated. 
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Algorithm  GRAPH 

Input:  An  rXc  totally  nonseparable  {0,l}-matrix  M .  We  assume  (see  the  first  paragraph  of  sec¬ 
tion  4)  that  the  first  column  of  M  is  a  singleton.  (Such  a  column  can  always  be  added  without 
affecting  the  realizability  of  M.) 

Output:  The  conclusion  that  M  is  not  realizable,  or  a  realizing  graph  G  for  M.  We  use  the  nota¬ 
tion  introduced  in  section  3. 

Step  Gl.  Mi  is  realized  by  a  bond  Gx  with  two  edges.  Set  {Gj}  where  for  pm(Gj)  we 
choose  the  non-tree  edge  of  Gi.  (Note  that  D\  is  a  decomposition,  but  not  a  t- 
decomposition  because  Gx  has  only  two  edges.  All  further  Dj,  j  >  2,  will  be  t- 
decompositions.)  If  c  =1,  stop  ( M  is  graphic).  Otherwise,  set  j  <—  2. 

Step  G2.  Set  P*~Pj  and  D*—Dj_v  Apply  HYPOPATH.  If  P  is  not  a  hypopath  of  m(D), 
stop. 

Comment.  If  a  maximal  graphic  subset  of  columns  is  desired,  not  just  a  determination  of 
the  realizability  of  M,  then  instead  of  terminating  when  P  is  not  a  hypopath,  simply  discard 
column  j  and  continue  the  algorithm. 

Step  G3.  Set  C*—Cj  and  apply  UPDATE.  (Note  that  the  quantities  Ky  K&  Uy  «2  required 
for  input  to  UPDATE  are  calculated  in  HYPOPATH,  and  the  associated  calls  to  TYPING, 
through  application  of  rules  (R1)-(R5).)  Set  Dj+—D‘. 

Step  G4-  If  j  =  c  ,  set  G  <—  m(D})  and  stop;  otherwise  set  and  go  to  Step  G2. 


The  validity  of  GRAPH  follows  from  the  results  of  sections  4  and  5. 


DATA  STRUCTURES.  Each  of  the  columns  and  rows  of  matrix  M  is  assumed  to  have  a  distinct 
name,  taken  from  the  integers  l,...,r  +  c;  all  marker  edges  created  during  the  algorithm  are  also 
given  distinct  names,  one  for  each  of  the  two  graphs  containing  them.  Each  member  of  each 
decomposition  is  named,  and  each  node  of  each  member  is  named.  These  node  names  for  the 
graphs  in  any  fixed  decomposition  will  always  be  distinct,  even  though  after  merging  several  may 
be  identified. 

(Dl)  M:  The  given  matrix  is  assumed  stored  in  column  list  form.  Thus,  for  each  column  we  can 
access  in  linear  time  the  list  of  rows  with  ones  in  this  column.  For  each  row  name  we  keep 
an  indicator  to  designate  if  this  row  is  in  a  previously  processed  column. 

(D2)  Members  of  D,  where  D  is  the  current  t-decomposition:  These  are  stored  and  accessed  using 
the  standard  disjoint  set  union  data  structure  with  path  compression  (see  Tarjan  [15]). 
Thus,  for  given  distinct  member  names  x,y  we  can  perform  makeset(x)  (assuming  x  is  not 
yet  in  any  set,  create  a  single  element  set  {ar}  with  name  ar),  }ind[x )  (find  the  name  of  the 
set  containing  x)  and  link(x,y)  (form  the  union  of  the  two  sets  containing  x  and  y).  The 
sets  are  stored  as  arborescences,  the  root  of  which  is  the  name  of  that  set.  Computing 
find(x)  involves  tracing  a  directed  path  from  x  to  the  root  of  its  tree,  and  then  “compress¬ 
ing”  this  path  so  that  all  nodes  on  it  point  directly  to  the  root.  The  operation  link{x,  y )  is 
performed  by  first  performing  find(x)  and  find(y),  and  completed  by  creating  an  arc 
pointing  from  the  root  of  the  shallower  of  the  two  trees  to  that  of  the  deeper. 
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Nodes  of  members  of  D:  Disjoint  set  data  structure,  as  in  (D2). 

Edge  locations :  Pointers  from  edges  to  the  names  of  the  members  containing  them.  Thus, 
to  find  the  member  containing  an  edge,  we  first  access  the  member  x  pointed  to  by  the 
edge,  and  then  perform  a  (D2)  find(x). 

Predecessor  function  p  of  D :  Pointers  from  member  names  to  member  names. 

Member  designation :  For  each  member  name,  a  indicator  designating  whether  it  is  a  bond, 
polygon,  or  prime. 

Polygon  sizes :  The  number  of  edges  in  each  polygon. 

Polygon  edge-sets:  Doubly-linked  lists. 

Edge  ends:  Pointers  from  each  edge  to  its  ends.  Note  that  for  a  given  edge,  once  the  two 
ends  u ,  v  have  been  accessed  it  will  still  be  necessary  to  perform  two  (D3)  finds  in  order  to 
find  the  current  names  of  these  ends.  For  marker  edges,  one  of  the  two  ends  is  designated  + 
and  one  -  .  This  serves  to  orient  D . 

(DIO) Parent  markers:  pm(x)  for  each  member  name  x. 

(Dll)  Child  markers:  For  each  current  member  x  ■=£  root(D),  pm(x)  is  associated  with  some  child 
marker  of  find(p(x)).  We  keep  a  pointer  from  x  to  this  child  marker.  In  the  case  when  x  is 
a  bond  and  find(p(x))  is  a  polygon,  we  also  keep  a  pointer  from  the  child  marker  to  x. 

DATA  STRUCTURE  APPLICATION  AND  UPDATE.  We  discuss  here  the  highlights  of  how 

(Dl)-(Dll)  are  applied  and  updated. 

(6.1)  Executing  Step  Hi:  P  and  C-P  are  computed  in  time  0(\C\)  using  (Dl),  and  the  indica¬ 
tors  for  C-P  are  updated.  Then  using  (D4)  and  (D2)  we  construct  a  list  of  members  in  D 
containing  edges  of  P.  This  requires  0(|P|)  (D2)  finds.  Using  this  list  and  (D5),  the  com¬ 
putation  of  D  and  simultaneously  the  depth  partition  (jt0,  ...,  tt,)  required  as  input  by  TYP¬ 
ING  uses  0(\D  |)  additional  (D2)  finds. 

(6.2)  Executing  Step  Tl:  If  H  is  a  bond,  as  determined  by  (D6),  this  is  straightforward.  Other¬ 
wise,  let  X  =  P  (^)E(H).  A  list  of  the  elements  in  X  can  be  conveniently  created  for  each 
HeD  as  D  is  computed.  If  H  is  prime,  then  using  (D9)  and  (D3)  we  construct  a  representa¬ 
tion  for  H[X]  that  uses  space  0(  \  X  | ),  and  do  this  with  0(  \  X  | )  (D3)  finds. 

The  most  complicated  case  occurs  when  H  is  a  polygon  and  RELINK1  is  applied.  Start¬ 
ing  with  the  marker  m  =  pm(H),  we  swap  edges  in  pairs  in  order  to  build  a  path  L  in  H 
one  end  of  which  is  m,  and  every  other  of  which  is  in  X.  At  a  general  step  we  take  the 


(D3) 

(D4) 

(D5) 

(D6) 

(D7) 

(D8) 

(D9) 
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non-m  end-edge  of  L  and  find  its  neighbor.  If  this  edge  is  not  in  X,  it  is  swapped  in  H  with 
the  next  edge  of  X  not  already  in  L.  This  uses  (D3),  (D8)  and  (D9)  and  updates  (D8)  and 
(D9).  It  requires  time  0(  |X| )  (D3)  finds. 

(6.3)  Executing  Steps  T3  and  H3:  In  addition  to  the  work  in  (6.2),  (D9)  is  applied  to  accomplish 
the  reorientation.  See  (4.7). 

(6.4)  Executing  (R1)-(R5):  The  application  of  (R5)  requires  special  consideration  since  the  final 
assignment  statement  may  require  accessing  a  child  of  a  node.  However,  this  happens  only 
for  polygon-parent,  bond-child  pairs  and  (Dll)  may  be  used. 

(6.5)  Executing  Step  Ul:  Determining  if  Kx  is  a  polygon  takes  constant  time  because  of  (D6). 
(Ul.l)  takes  constant  time,  as  does  (U1.2).  In  (Ul.3)  particular  care  must  be  taken  with 
(D4).  Thus,  when  two  polygons  are  formed  from  one,  then  only  one  of  these  can  retain  the 
name  of  the  old  polygon  member,  and  hence  for  every  edge  in  one  of  the  new  polygons  (D4) 
must  be  updated.  In  this  update  we  always  choose  the  piece  Ly  that  meets  the  set 
W  =  (P  p]  E(KX))  U  F  where  F  is  the  set  of  child  markers  of  Kx  such  that  P  meets  the 
corresponding  complete  child.  It  is  easy  to  see  then  that  Ly  C  W,  and  hence  that  0(|P|) 
work  will  accomplish  the  update. 

We  must  also  update  (D5)  for  all  children  of  the  polygon  Ly(_J{/y}  having  the  same 
name.  This  can  be  done  in  0(\D  |)  time  since  Ly  C  W  implies  all  children  are  in  D. 

(6.6)  Executing  Step  U2:  First  R  must  be  calculated,  which  uses  0(|i?  |)  (D2)  finds.  (These  will 
be  constant-time  finds  given  that  the  reduced  ^decomposition  for  P  has  already  been  found 
in  HYPOPATH.)  (U2.1)  requires  four  (D3)  finds  to  check  adjacency  of  my  and  my_j.  For 
(U2.2)  and  (U2.3),  the  nontrivial  work  involves  the  applications  of  SQUEEZE.  Again,  by 
the  convention  introduced  in  (6.5),  the  total  work  for  the  SQUEEZES  involves  0(|P|) 
constant-time  updates  in  (D4),  and  0(\R  ])  constant-time  operations  on  each  of  (D2),  (D3), 
(D5)-(D11).  Finally,  in  (U2.4)  the  forming  of  m(R)  requires  0(\R  |)  (D2)  and  (D3)  links. 

(6.7)  We  remark  that  (D7)  is  not  necessary  to  implement  the  algorithm.  However,  its  availability 
will  sometimes  speed  up  the  computation  of  T(H)  when  H  is  a  polygon. 
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7.  TIME  AND  SPACE  BOUNDS 

We  begin  by  giving  bounds  for  one  application  of  HYPOPATH  and  one  application  of 
UPDATE.  Note  that  (Rl)-(R5)  of  section  5  are  applied  in  HYPOPATH. 

(7.1)  Lemma.  One  application  of  HYPOPATH  uses  time  bounded  by  the  time  required  for 
0(\P\  +  |  D  | )  find  operations  using  (D2)  and  (D3). 

Proof.  Step  HI  takes  0(|P|  +  |D|)  (D2)  finds  by  (6.1).  Now  consider  Steps  Tl,  T3  and  H3. 
For  QGD,  let  be  a  list  of  the  complete  children  of  Q  in  D  (this  list  may  be  empty). 

Define  Wq  —  (P  p|  £(<9))  LJ  {wh,...,mt}  where  m,- =  pm(H{)  (1  <  i  <  t),  and  set 
H—Q/[Hi,...,Ht\.  By  (6.2)  and  (6.3),  T(H)  and  an  orientation  f'  such  that  is  good 

can  be  found  with  0(  |  Wq  | )  (D3)  finds.  But 

£  I  I  =  |F|  +  |D|-1. 

Qet> 

This  proves  the  result. 

(7.2)  Lemma.  One  application  of  UPDATE  uses  time  bounded  by  0(  |  G\  +  |i?  |)  finds  using 
(D2)  and  (D3)  (R  is  defined  in  Step  U2). 

Proof.  Step  UO  can  clearly  be  carried  out  in  time  0{  \  C  | ).  The  remainder  of  the  proof  follows 
from  (6.5)  and  (6.6). 


For  the  final  derivation  of  the  time  and  space  bounds  we  need  the  results  (7.3)-(7.1l). 

(7.3)  Lemma.  Let  G  be  a  nonseparable  graph  with  r  >3  nodes.  Let  D  be  a  (-decomposition  of 
G  and  let  be  the  prime  members  of  D.  Then 

£|W(//,-)  I  <  2r-4. 

1=1 


Proof.  The  result  follows  easily  by  induction  on  r. 
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The  algorithm  GRAPH  proceeds  by  producing  a  sequence  of  ^decompositions.  Let 
Di,...,Dc  be  the  sequence  produced  for  a  given  rXc  {0, l}-matrix.  For  each  index  j  (l</<  c) 
such  that  Step  U2  is  applied,  denote  by  Rj  the  path  R  of  nodes  merged  in  (U2.4).  Every  prime 
node  H  that  occurs  in  the  algorithm  originates  from  the  application  of  (U2.4)  to  some  Rj)  H  may 
then  be  further  modified  by  application  of  (Ul.l)  or  (U2.1),  and  so  several  H s  may  have  the  same 
associated  j ,  but  the  value  of  j  for  a  given  H  is  still  well  defined.  Now  define,  recursively,  the 

sets  IH  by  IH  =  {;}  (J  /  (J  l\. 

I  J  prime  j 

(7.4)  Lemma.  For  each  prime  H  that  occurs  as  a  node  of  some  decomposition  in  GRAPH, 

\N(H)\  >  3  +  £(|R,.|-F|R,|/2l). 

iei„ 

(M  is  the  least  integer  not  less  than  x.) 

Proof.  Letp,  =  |R,  |  (*€///).  The  proof  is  by  induction  in  If  Ijj  —  {i },  then 

\N(H)\  =  2(1 -pf)  +  E  \  W)  I- 

If  /  is  a  bond,  then  |  N(J)  \  —  2,  and  if  J  is  not  a  bond  then  |  N(J)  \  >  3.  Since  no  two  bonds  can 
share  a  marker  edge,  and  neither  end-graph  of  Rj  can  be  a  bond,  by  Lemma  5.1,  the  number  of 
bonds  in  Rj  is  bounded  by  fp,  /2]  -1.  Therefore 

|  N(H)  |  >  2(l-pi)  +  3 p,.  -  r  Pi/21  +  1=3  +  p,-  -  \Pi/2\  . 

Now  let  |  Iff  j  1.  Define  tn  =  max{  1 1  ^  Iff }  and  Is  ==  In-{m}.  For  each  prime  node  J£.Rm, 
induction  yields 

\N(J)\  >  3  +  E(p.-rp./2i). 

ieij 

Let  Rm'  =  {J  ERm.  J  prime}  and  Rm"  =Rm—Rm'.  Then 

\N(H)\  >  2(1  -pm)+  E  (3+E(P.-rP./2l)l  +  E  1^)1- 

/efim(  iei,  )  J€Rm- 

Since  {Ij\  J £Rm'  }  partitions  IH' ,  it  follows  that 
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\mm  >  2(i-Pm)  +  3|/?m  i+  s(p,-fp,/2i)  +  3(pm-|i?m'i)-(rpm/2i-i) 

•w 

= 3  +  s  ?.■  -  s  rp,/2i 

•€/h  >eiH 


(7.5)  Theorem.  For  r  >  3, 

E  Kl  <  6r-12- 

aii  /?, 

Proof.  Let  P  be  the  set  of  prime  members  of  Dc.  By  Lemmas  7.3  and  7.4, 

6r-12>  3£(3+E(|i?,|-r|i?,|/2l) 

H€p{  ieiH 

>  E  (9  +  E  i^.l)  >  E  Kl- 

HeP  (  i£lH  )  allRt 

Theorem  7.5  says  that  the  total  number  of  graphs  ever  involved  in  a  (U2.4)  merge  is  0{r). 
For  each  t-decomposition  D{,  let  D(  be  the  corresponding  reduced  ^decomposition.  Then 
RiC.  D{. 

Next  we  bound  the  number  of  nodes  in  Dit  but  not  in 

(7.6)  Lemma.  Let  H  be  a  nonseparable  graph  and  let  X  be  a  cycle  of  H.  Let  D  be  a  t- 
decomposition  of  H  and  let  D  be  the  reduced  (-decomposition  with  respect  to  X.  Then 
\D\  <  |X|+1. 

Proof.  The  theorem  is  clear  if  |  D  \  —  1.  If  D  has  depth  one,  then  all  but  at  most  one  node  of  D 
contains  edges  of  X.  Again  the  result  is  clear.  Hence,  we  may  assume  that  there  is  a  pendant 
node  K  of  D  where  K  is  at  least  distance  two  from  root(D). 

Let  K'  be  the  parent  of  K,  and  let  X'  =(X-E{K)){J  {pm(K)}.  Clearly  X'  is  a  cycle  of 
m(D'),  where  D'  is  D  with  K  deleted.  If  K'  has  more  than  one  child,  D'  is  a  reduced  t- 
decomposition  for  X' .  In  this  case,  the  theorem  follows  by  induction  on  \D\-  Suppose  K'  has 
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only  K  as  a  child  and  that  Xf") E{K')  =  <j>  (otherwise,  again  we  may  apply  induction).  But  then 
pm(K)  is  parallel  to  pm(K')  in  K'  and  so  K'  is  a  bond. 

Now  by  the  choice  of  K,  K'  is  not  root(D);  moreover,  p{K')  is  not  a  bond.  Hence  D"  ,D 
with  K  and  K'  deleted,  is  a  reduced  (-decomposition  for  X' .  In  addition,  note  that 
|Xp)£(X)|  >2;  otherwise  K  has  and  edge  parallel  to  its  parent  marker,  and  is  not  a  bond. 
Since  at  least  two  edges  of  X  have  been  deleted,  and  exactly  two  nodes  of  D,  the  result  follows  by 
induction. 

(7.7)  Theorem.  Let  DU...,DC  be  the  sequence  of  reduced  (-decompositions  produced  by  GRAPH. 
Then 

Y  |  A  I  <  2n  +  c+6r-12. 

i=i 

Proof.  Consider  a  particular  D ,.  Let  /?,•  —  {Ki}  if  Kl  =  K2  in  UPDATE.  Now  D{  can  be 
viewed  as  a  collection  of  trees  attached  to  R{.  For  a  particular  one  of  these  trees  T' ,  let  rn'  be 
the  marker  linking  T'  to  Then  clearly  X=(P,-+iP)£l(m(T''))U  {m'}  is  a  cycle  in  m(T'), 
and,  indeed,  T'  is  a  reduced  (-decomposition  for  X.  Hence,  by  Lemma  7.6, 
\T'  |  <  |X|+1<2|X|.  Summing  over  all  such  T' ,  we  have  that  the  to#al  number  of  nodes  of 
Di  not  in  R{  is  bounded  by  2  |P1+1 1 .  Hence, 

E  I A  I  <  E(l^l+2|Pt+1|)  <  c+6r-12  +  2n 

i=i  i=i 

by  Theorem  7.5. 

(7.8)  Lemma.  Y  I  Cy— Py  |  <  2r  . 

\C]-PJ\>2 

(7.9)  Lemma.  For  a  single  application  of  GRAPH,  the  total  number  of  node  names  created  for 
nodes  of  members  of  (-decompositions  and  the  total  number  of  member  names  created  are  both 
bounded  by  0{r). 

Proof.  For  a  (-decomposition  D,  define  S(D)=  Y  ( I  C  |  -3:C  is  the  edge-set  of  a  polygon  of  D). 
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Note  that  in  order  for  (U1.3)  or  SQUEEZE  to  be  applied  it  must  be  that  S(D)  >  0.  Moreover,  the 
application  of  either  of  these  operations  reduces  S{D).  On  the  other  hand,  S(D)  can  only  be 
increased  in  Step  U0,  and  then  only  by  |  G-P  |  -2,  whenever  |  G—P  |  >  2.  Hence,  it  follows  from 
Lemma  7.8  that  (Ul.3)  and  the  SQUEEZE  procedure  can  be  applied  at  most  O(r)  times. 

Now  define  S'(D)  =  ^  ( |  C  |  - k(C):  G  is  the  edge-set  of  a  polygon  of  D)  where  k(C)  is  the 
number  of  markers  in  G  that  are  also  in  some  bond.  Note  that  in  order  for  (Ul.2)  to  be  applied  it 
must  be  that  S'  (D)  >  0.  S'(D)  can  only  be  increased  by  application  of  Step  U0  or  the  SQUEEZE 
procedure.  The  total  increase  over  all  applications  of  Step  U0  is  O(r)  by  Lemma  7.8.  On  the 
other  hand,  each  application  of  SQUEEZE  increases  S'(D)  by  only  a  small  constant.  In  view  of 
the  calculations  of  the  previous  paragraph,  the  total  increase  over  one  application  of  GRAPH  is 
therefore  bounded  by  O(r).  But  each  application  of  (Ul.2)  decreases  S' (D).  It  follows  that 
(Ul.2)  can  be  applied  at  most  O(r)  times. 

By  the  calculations  of  the  previous  two  paragraphs  we  may  ignore  (Ul.2),  (Ul.3),  (U2.2)  and 
(U2.3)  in  counting  the  number  of  nodes  and  the  number  of  members  created.  The  number  of 
applications  of  (U2.1)  is  bounded  by  0(r)  because  of  Theorem  7.5.  (U2.4)  adds  no  node  or 
member  names  (the  new  member  m(R)  inherits  its  name,  through  (D2)  links,  from  one  of  the 
members  in  R).  (Ul.l)  also  adds  no  node  or  member  names.  Finally,  the  number  of  node  and 
member  names  created  in  Step  U0  is  bounded  by  O(r)  by  Lemma  7.8. 

(7.10)  Lemma.  The  total  number  of  edge  names  created  by  an  application  of  GRAPH  is 
bounded  by  0(r  +  c). 

Proof.  A  direct  examination  of  UPDATE  shows  that  no  more  than  |(7-P|+6  edges  can  be 
added  by  any  one  application. 

• 

We  now  introduce  a  result  of  Tarjan  [15].  Define  the  so-called  Ackerman  function  A(i,j)  on 
pairs  of  positive  integers  i,j  by:  A(l,/)  =  2;,  A(t',l)  =  A(i - 1,2)  for  »>  2  and 

-A(*,y)  =  A(i-l,A{i,j -\))  for  *,/ >  2.  For  integers  m>n>  1,  define 
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a(m,  n)  =  min{«  >  1:J4(«,  [m/nj  )  >log2n}.  The  function  a(m,n )  is  very  slow  growing,  being  for 
all  practical  purposes  never  bigger  than  4.  (>4(4, l)  =  /(17),  where  /( 1)  =  2  and  f(k)  =  2^k~1'1  for 
k  >2.  Thus, /(4)  =  266636.) 

(7.11)  Theorem.  (Tarjan  [14])  The  disjoint  set  union  algorithm  with  path  compression  runs  in 
time  0(ma(m  ,n))  where  m  is  the  total  number  of  operations  (makesets,  links  and  finds)  and  n  is 
the  number  of  elements  in  the  underlying  set. 

(7.12)  Theorem.  Given  an  rXc  (0,  l}-matrix  with  n  nonzero  entries,  GRAPH  runs  in  time 
0(na(n,r))  and  uses  space  0(n). 

Proof.  The  time  bound  for  Steps  G1-G3  follows  from  (7.1),  (7.2),  (7.7),  (7.9)  and  (7.11).  For 
Step  G4  note  that  G  can  be  computed  directly  using  (D9)  and  0(c)  (D3)  finds.  The  space  bound 
follows  from  (7.9),  (7.10)  and  an  examination  of  (Dl)-(Dll).  The  determining  factor  is  the  storage 
for  the  input  matrix.  We  note  that  some  temporary  storage  is  involved  in  the  implementation  of 
several  of  the  steps,  principally  Steps  Tl,  T3,  and  H3  (see  (6.2)).  This  storage  is  not  recorded  in 
(Dl)-(Dll).  However,  it  can  easily  be  verified  never  to  exceed  0( \P  |  +  \D  | ). 

It  is  open  whether  there  is  an  0(n )  algorithm  for  graph  realization.  Recently  Gabow  and 
Tarjan  [8]  have  given  some  conditions  under  which  disjoint  sets  can  be  treated  in  linear  time,  and 
although  such  a  reduction  in  the  time  required  for  the  applicable  parts  of  GRAPH  would  imply  a 
linear  time  bound,  our  problem  does  not  seem  to  meet  their  conditions.  Indeed,  we  conjecture 
that  there  is  no  linear-time  (i.e.,  O (n )j  algorithm  for  graph  realization. 
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APPENDIX.  AN  EXAMPLE 

Consider  the  9X8  matrix  M  in  Figure  3: 

10  11  12  13  14  15  16  17 

1  11000001 
2  00111001 

3  00111001 

4  00111011 

5  0  0  0  0  0  1  1  1 

6  0  0  1  1  1  1  0  0 

7  0  1  1  1  1  0  0  0 

8  01101000 

9  01000000 

Figure  3:  Input  Matrix  M 


M  is  totally  nonseparable  since  each  column,  except  the  first,  has  a  1  in  a  row  with  a  1  in  a 
previous  column.  We  use  the  labels  BltB2,  ...  for  the  graphs  that  occur  in  t-decompositions  dur¬ 
ing  the  execution  of  the  algorithm. 


Column  11:  In  Step  G2  we  have  P  =  {1}.  Note  that  D=Dl  is  not  a  t-decomposition  since  B j 
has  only  two  edges,  but  is  a  t-decomposition  in  all  other  respects. 

Step  HI  is  trivial  since  D  has  only  one  member.  In  Step  H3,  t  =  0,  so  H—  Q  =Bl.  Bx  is 
by  convention  not  a  polygon,  since  it  has  two  edges,  and  so  RELINK2(i?i)  is  not  applied.  Note 
that  |  P  |  =1,  and  so  P  is  automatically  a  path.  Note  also  that  (R4)  and  (R5)  are  applied  in  Step 
H3.  Since  TYPING  was  not  used  in  this  case,  (Rl)-(R3)  were  not  applied.  Thus,  K1=K2  =  Bx. 
No  change  occurs  in  (R5)  since  Bx  is  a  bond,  even  though  Ki=K2  and  the  ends  of  P  are  the  ends 
of  =  10. 
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In  Step  G3  we  have  C=  {1,7,8,9,11}-  The  identities  of  u1,u2  as  required  by  the  input  for 
UPDATE  are  clear  since  Bx  has  only  two  nodes.  In  Step  UO  we  have  G-P  =  {7,8,9,11}-  Let 
/=18.  Adding  a  polygon  B2  with  edges  {7,8,9,11,18}  to  D'  =  {B1}  we  obtain  the  non¬ 
decomposition  D'  as  in  Figure  5  (“non”  because  there  is  no  marker  edge  shared  by  Bx  and  B2). 

10 

Bi  <^> 

18 


Figure  5s  £>*  for  column  11 


In  Step  Ul  we  apply  (Ul.l)  since,  by  convention,  Bx  is  not  a  polygon.  This  yields  the  t- 
decomposition  D2  in  Figure  6.  Note  that  the  marker  edge  18  has  been  oriented.  The  orientation 
is  arbitrary,  as  suggested  by  the  discussion  following  the  proof  of  Lemma  5.1.  Note  also  that  edge 
18  is  given  the  same  name  in  both  Bx  and  B2.  This  causes  no  difficulty  here,  in  contrast  to  the 
situation  in  a  computer  implementation,  as  the  distinct  identities  of  these  two  edges  is  clear  from 
the  picture. 

B  A' 

Bi  » 

sfr' 

j 

B2  11 


18  U1 
O’ 


Figure  8:  D2  (column  11) 


Column  12:  We  now  have  P={7,8}.  In  Step  HI  the  reduced  decomposition  is  D  —  {B2}  since 
this  member  of  D  contains  all  the  edges  of  P.  In  Step  H3,  f=0,  so  H  =  Q  —  B2.  Since  B2  is  a 
polygon  we  apply  RELINK2(B2)-  Now  Wq  =P  =  {7,8}  (see  definition  of  Wq  preceding  Lemma 
4.3),  and  Z  =  <f>  since  t=  0  (see  definition  of  Z  preceding  RELINK1).  But  Wq-Z  =  { 7,8}  is 
already  a  path  in  B2,  and  so  nothing  happens.  Also  no  reorientation  occurs  in  Step  H3  because 
<  =  0,  or,  more  generally,  because  Z  =  <f>.  Because  of  (R4),  Kx  —  K2=  B2.  No  change  occurs  in 
(R5)  since  and  u2  are  not  adjacent  (see  Figure  6). 
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In  Step  G3  we  have  <7=  {2,3,4,6,7,8,12},  and  in  Step  UO  C-P  =  {2,3,4,6,12}.  Let  /  =  19. 
Adding  a  polygon  63  with  edges  {2,3,4,6,12,19}  to  D'  we  obtain  the  non-decomposition  D*  as  in 
Figure  7. 


Figure  7:  D*  for  column  12 


In  Step  U1  we  apply  (U1.3).  Letting  / j  =  20  and  / 2=  21,  and  where  the  path  Ly  has  edges 
{9,11,18}  and  62  has  edges  {7,8},  we  obtain  the  {-decomposition  Z)3  in  Figure  8. 
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Figure  8:  D3  (column  12) 

Column  13:  We  now  have  P  =  {2, 3, 4, 6, 7}.  In  Step  HI  the  reduced  decomposition  is 
D  =  {.63,54,135},  as  given  in  Figure  9,  where  the  inclusion  of  64  is  forced  by  connection.  This 
leads  to  our  first  application  of  TYPING. 

y’ 


Figure  9:  D  for  column  13 
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The  depth  partition  is  given  by  tt0={.B4}  and  tti  =  {Bb,Bb},  so  that  s  =  1.  In  Step  Tl  we 
must  apply  both  RELINKl(B3)  and  RELINK1(R6)  since  Bz,Bb  are  both  polygons.  However,  in 
both  cases  Z  =  <f>\  moreover,  WBjs  =  {2, 3, 4, 6}  and  WB&={7}  are  both  paths  incident  to  the  parent 
markers,  19  and  21,  respectively.  Hence,  no  relinking  actually  occurs,  and  we  have 
T(Bs)=T(Bb)  =  3. 


We  must  now  apply  Step  H3.  We  have  Q=Bi,  HX  =  BS  and  H2  —  B 6.  The  current  orien¬ 
tation  is  given  by  the  directions  on  the  marker  edges  19  and  21.  Since  Q  is  not  a  polygon  we  do 
not  apply  RELINK2.  It  remains  only  to  determine  a  new  orientation  /.  As  suggested  in  (4.7),  to 
do  this  all  we  need  to  know  about  and  H2  is  that  both  are  type  3,  in  addition  to  knowing 
which  parent  marker  end-nodes  meet  WH i  and  WH^,  respectively.  These  are  the  nodes  x  and  y  in 
Figure  9.  Since  their  corresponding  nodes  on  Q  under  /,  x'  and  y’ ,  are  not  equal,  one  of  19,21 
must  be  reoriented.  We  reorient  19,  which  yields  the  reduced  decomposition  D  given  in  Figure 
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Figure  10:  Reoriented  D  for  column  13 


In  Step  G3  we  have  G—  {2,3,4,6,7,13}.  For  the  input  requirements  of  UPDATE  we  note 
that  Ky,K2 ,«i,u2  are  chosen  in  Step  Tl  by  the  application  of  rule  (Rl).  Thus,  we  have  Kl  =  Bb 
and  K2  =  Bb  (K2  =  B3,  Ki  =  Bb  is  an  equally  valid  choice),  where  u ltu2  are  specified  in  Figure 
10.  In  Step  U0  we  have  D"  as  in  Figure  8,  except  that  edge  19  has  been  reoriented  in  P4.  Since 
C-P  =  {13}  we  obtain  /  =  13.  Nothing  is  added  to  D’  in  this  step. 

Now  we  come  to  Step  U2.  For  the  path  R  we  have  J1  =  Bb,  J2  =  Bit  Jz  —  Bb,  s—  2, 
m!=19,  and  m2  =  21.  Note  that  S4  is  a  bond,  and  that  p(R4)  =  P2  is  n°l  in  B.  However,  J94  has 
only  three  edges.  Thus,  no  change  occurs  in  (U2.1).  In  (U.2)  we  apply  SQUEEZE(L)  for 
L  ==  (2, 3, 4, 6};  all  other  Ls  encountered  have  only  one  edge  so  that  SQUEEZE  has  no  effect  in 
these  cases.  The  resulting  D*  is  given  in  Figure  11. 
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Figure  11:  R  for  column  13 


Nothing  occurs  in  (U2.3)  since  R  has  no  internal  polygons.  Finally,  applying  (U2.4)  we 
obtain  the  i-decomposition  Z)4  of  Figure  12,  where  B&  is  the  graph  G  produced  by  the  merge. 
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Figure  12:  Z)4  (column  13) 


Column  14:  We  have  P  =  (2, 3, 4, 6, 7, 8},  D  =  {BS,B6}  7r0={P3},  and  n x  —  {Be}.  In  Step  Tl, 
RELINK(P6)  has  no  effect,  and  T(Be)—  1.  This  brings  us  to  Step  H3  with  Q  =B3,  HX=B6  and 
t  =  1.  Since  T(HX)  =  1  no  reorientation  occurs.  KX  =  K2  =  B3  because  of  rule  (R4)  applied  in 
Step  H3,  with  ux,u2  as  indicated  in  Figure  12.  In  Step  UO,  /  =  14.  Finally,  (Ul.l)  is  applied 
yielding  Db  as  in  Figure  13. 


Column  15:  We  have  P  =  {6}  and  D  =  {Pg}-  Step  H3,  RELINK2(P6)  switches  the  positions  of 
4  and  6,  where  Wq  =  { 6},  Z  =  <j>,  so  that  6  is  adjacent  to  pm(Be)  —  22.  Because  <=0,  no 
reorientation  occurs,  KX  =  K2=B6,  by  (R4),  and  ux,u2  are  the  ends  of  6  in  P6.  In  Step  UO  we 
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take  /  =23  (a  new  element)  and  introduce  a  polygon  B 7  on  edge-set  {/ }  {_J  (G-P)  =  {5,15,23}. 
Next  (U1.2)  is  applied,  with  /'  =6.  Let  f"  =24.  Replacing  edge  6  by  edge  24  in  Be,  and  adding 
a  bond  Bs  on  edges  {6,23,24}  we  obtain  De  as  in  Figure  14. 


b6  b8 


Figure  14s  D6  (column  15) 


Column  16:  We  have  P  =  {4,5},  D  =  {Bft,B7,Bs},  7r0={B6},  Tri  —  {BS},  and  7r2={B7}.  Applying 
Step  T1  to  B7  we  see  that  RELINK1(B7)  has  no  effect,  and  determine  that  T(J37)  =  3,  Kl—B7 
(by  rule  (Rl)),  and  that  is  the  node  of  B7  not  incident  to  edge  23.  Next  we  apply  Step  T3  with 
Q=BS  and  Hl  =  B7  (t  =  1).  No  reorientation  occurs,  and  we  deduce  that  T(Q[Hi})  =  3.  This 
completes  the  application  of  TYPING,  and  brings  us  to  Step  H3.  In  applying  RELINK2(B6),  we 
note  that  Wq  =  {4,24}  and  Z  =  { 24},  and  4  is  switched  with  3  to  make  it  adjacent  to  24.  The 
application  of  rule  (R4)  yields  K2  —  B6  and  u2  equal  to  the  end-node  of  4  not  shared  by  24. 

Now  in  Step  UO  we  have  C—  {4,5,16}  and  /  =  16.  (U2.1)  has  no  effect.  (U2.2)  results  in 
squeezing  the  path  {2,3,22}  from  B$.  (U2.3)  has  no  effect  because  R  has  no  internal  polygons. 
Finally,  applying  (U2.4)  we  obtain  D7  as  in  Figure  15. 


9 


37 


Column  17:  We  have  P  =  {1,2, 3, 4, 5},  D  =  D7,  and  blocks  in  the  depth  partition 
7r0=  (Bj },..., tt4=  {Bg}.  In  TYPING  we  obtain  T{H)  =  3  for  each  application  of  Steps  T1  and 
T3.  RELINK  1  is  applied  twice,  once  to  Be  and  once  to  B2,  but  in  neither  case  results  in  any 
change.  The  first  application  of  Step  Tl  results  in  identifying  Ki  as  B%  and  as  the  node  of  Bs 
common  to  edges  5  and  15.  K2  is  chosen  to  B2  in  Step  H3  since  u2  is  common  to  edges  11  and 
18.  No  reorientations  occur  in  TYPING  or  HYPOPATH.  Finally,  in  UPDATE  Step  U2  is 
applied.  (U2.2)  results  in  squeezing  the  path  {11,9}  from  B2,  and  in  (U2.3)  the  path  {2,3}  is 
squeezed.  The  application  of  (U2.4)  results  in  Ds  as  in  Figure  16.  m(D8),  a  graph  realizing  the 
original  M,  is  given  in  Figure  17. 


Figure  16:  £>g  (column  17) 
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Figure  17:  A  realization  for  M 
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